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By using bulk magnetization, electron spin resonance (ESR), heat capacity, and neutron scattering 
techniques, we characterize the thermodynamic and quantum phase diagrams of Ba3Cr2 08. Our 
ESR measurements indicate that the low field paramagnetic ground state is a mixed state of the 
singlet and the 8^=0 triplet for H J_ c. This suggests the presence of an intra-dimer Dzyaloshinsky- 
Moriya (DM) interaction with a DM vector perpendicular to the c-axis. 

PACS numbers: 75.10.Jm, 75.25.+Z 



Over the last two decades, several exotic collective phe- 
nomena have been found in quantum (s = 1/2) mag- 
nets [H- One of the well-known examples is the con- 
densation of magnons induced by application of an ex- 
ternal magnetic field [H, [1]. The field-induced quan- 
tum phase transition has been experimentally observed in 
coupled quantum dimer systems, such as TlCuCla [1, 
and BaCuSi206 [5]. If the field is applied along an axis 
of U(l) symmetry, the induced quantum critical point 
(QCP) belongs to the Bose-Einstein condensation (BEC) 
universality class [6]. One signature of the BEC univer- 
sality class is the power law for the increase in the transi- 
tion temperature, Tc close enough to the critical field Hd 
: Tc (X [H — HciY^^ ^ with the critical exponent (j) = d/2 
for spatial dimensionality d > 2 |6|]. The BEC value 
of (j) has been found in the two-dimensionally (2d) cou- 
pled dimer system BaCuSi206 [s^, a frustrated 2d antifer- 
romagnet CS2CUCI4 [7], and a three-dimensionally cou- 
pled s = 1 chain NiCl2-4SC(NH2)2 0. However, values 
of <j) higher than 3/2 have been reported for the three- 
dimensionally coupled quantum dimer compounds such 
as TlCuCls (0 = 2.2(1)) and KCuCls {(j) = 2.3(1)) for 
1 K < T < 5 K. Quantum Monte Carlo (QMC) simula- 
tions of U(l) invariant models [1, have shown that (j) 
can be easily overestimated by fitting the non-universal 
regime. On the other hand, when a narrow enough tem- 
perature range AT is considered (AT ^ {).ATmax for an 
s = 1/2 XXZ magnet), the QMC results yield the ex- 
pected (j) = d/2 exponent [lO^. AT is the typical size of 
the quantum critical region. 

A question that arises is how robust is the BEC 
quantum critical behavior against the presence of mag- 
netic anisotropy. A case in point here is the three- 
dimensionally coupled quantum dimer system Ba3Cr204 
where the tetrahedrally co-ordinated Cr^+(3(i^;s = 1/2) 
form dimers along the crystallographic c-axis which in 



turn are coupled into a triangular lattice in the a6-plane. 
Unlike the related material, Ba3Mn208 that has the or- 
bitally nondegenerate Mn^+(3(i^;5 = 1) ions, Ba3Cr204 
undergoes a structural phase transition at 70 K due to 
the Cg Jahn- Teller active Cr^+ ion that leads to spatially 
anisotropic three-dimensional exchange interactions and 
lifts the frustration of the triangular lattice [ll| . Without 
an external magnetic field, Ba3Cr204 does not exhibit 
any long range order down to 1.5 K. When an external 
magnetic field is applied, however, the bulk magnetiza- 
tion data shows a field-induced transition into a magnet- 
ically ordered state at Hd 12 T and another transition 
into a fully polarized state at Hc2 ~ 23 T [l2| . 

In this letter, we report our systematic studies of 
the field-induced phase transitions in Ba3Cr208, using 
bulk magnetization, specific heat, electron spin resonance 
(ESR), and both elastic and inelastic neutron scattering 
techniques. Our ESR data showed that Ba3Cr208 re- 
alizes a system of three-dimensionally coupled quantum 
dimers with weak inplane Dzyaloshinsky-Moriya (DM) 
interactions less than 0.1 meV. Our specific heat and elas- 
tic neutron scattering data show that the H — T phase 
boundary of the field-induced canted antiferromagnetic 
(AFM) state near the lower critical field follows a power 
law Tc (x| Hc{T) - Hc{0) |^/^ that is consistent with 

= d/2 for 30mK < T < 1 K. 

Single crystals of Ba3Cr208 were grown at the Insti- 
tute for Solid State Physics (ISSP), University of Tokyo 
as well as at the University of Virginia, using fioat- 
ing zone image furnaces. Bulk magnetization, electron 
spin resonance and specific heat measurements were per- 
formed at the Institute for Materials Research (IMR) of 
Tohoku University using small single crystals, a 15 mg 
piece for H ± c and two co-mounted crystals of total 
weight of 20 mg for i^" || c. Elastic and inelastic neutron 
scattering measurements were carried out at the cold- 
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FIG. 1: (color online) (a) Bulk magnetization (M) as a func- 
tion of an external pulsed magnetic field (H) applied parallel 
and perpendicular to the c-axis, taken at T = 0.4 K. Ar- 
rows indicate upper and lower critical fields determined by 
the positions of the peaks that appear in the second deriva- 
tive, dPM/dH^^ as a function of H as shown in the inset, 
(d) Temperature dependence of the specific heat under var- 
ious constant steady magnetic fields along the c-axis, from 
H — 12.5 T to 17 T. The data are consistent with Aczel et 
al.^s recent results reported in Ref. p^ . 



neutron triple-axis spectrometer V2 at the Helmholtz- 
Zentrum Berlin (HZB). A larger single crystal of 2.5 g was 
used for the measurements and was oriented with (/i, I) 
as the scattering plane. A 15 T magnet with dilution in- 
sert was used. A flat analyzer configuration with a final 
neutron energy of Ef =3.7 and 5 meV and a horizon- 
tally focusing configuration with Ef = 5 meV were used 
for elastic and inelastic neutron scattering measurements, 
respectively. A cooled Be filter was placed before the an- 
alyzer to eliminate higher order neutron contaminations. 
Horizontal collimations were guide-60 '-open-open at V2. 

Figure [Tfa) shows bulk magnetization (My) measured 
using a pulsed magnet at T = 0.4 K and for fields up 
to 25 T. Upon ramping up, the ferromagnetic longitudi- 
nal component of the magnetic moment, My, increases 
rapidly at Hd = 12.15(5) T when ± c, as previously 
reported, and at Hd = 12.45(5) T when H \\ c. Upon 
further increasing My increases linearly up to a sat- 
urated value, Ms = gsjiB/Ci^^ at Hc2 = 23.05(5) T for 
± c and at Hc2 = 23.70(5) for H \\ c. The critical 
fields marked by the arrows are determined from the sec- 
ond derivative, cPM/dH^. Fig. [Tfb) shows specific heat 
data as a function of temperature measured at various 
fields up to 17 T applied along the c-axis. At 12.5 T just 
above i^ci, a weak A-type anomaly is observed at 0.94 K. 
As H increases up to 17 T, the anomaly becomes more 
prominent and the transition temperature also increases 
up to 2.84 K. 

In order to investigate the nature of the field-induced 
state above i^ci, we performed elastic neturon scatter- 
ing measurements as a function of magnetic field applied 
perpendicular to the c axis and temperature. As shown 
in Fig. [2](a), for H > Hd magnetic Bragg peaks ap- 
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FIG. 2: (color online) (a) Elastic neutron scattering data 
around the magnetic (0.5,0.5,3) reflection, measured at 0.15 
K and with three different magnetic fields, H = 9 T (black), 
12.5 T (blue), and 13.5 T (red circles). The inset shows the 
color contour map of the inelastic neutron scattering intensity 
of the Sz = 1 triplet excitation measured with H — 14.8 T. 

(b) The critical exponent 2/3 obtained from the data shown in 

(c) . See the text for detailed description, (c) if-dependence 
of the magnetic (0.5,0.5,3) Bragg reflection measured at var- 
ious temperatures from 30 mK to 2.16 K. H was perpendic- 
ular to the c-axis. The solid lines are fits to a power-law. 
(c) Schematic figure of the structure of the transverse mo- 
ments in the canted AFM phase of Ba3Cr208 for if ± c, 
obtained from our magnetic Bragg scattering data. Arrows 
represent the transverse magnetic moments along the c-axis. 
The magnetic unit cell is 2a x 2b x c. 



pear at Q = (0.5,0.5,/) with / = 3n (integer n). This, 
combined with the bulk magnetization data, clearly in- 
dicates that the field-induced phase is a canted antifer- 
romagnetic state. As H increases, the magnetic Bragg 
intensity increases. The if-dependence of the intensity 
has been measured at various Ts, as shown in Fig. ^c). 
Upon increasing T from 30 mK to 2.1 K, Hd increases 
from 11.95(2) T to 13.72(1) T. The if-dependent inten- 
sity is proportional to the square of the antiferromagnetic 
transverse component of the magnetic moment, M]_ and 
follows the power law, M]_ (x {H — HdY^ - The solid 
lines are the fits to the power law using all data avail- 
able, yielding 2/3 = 0.73(3) for 30 mK and 0.57(3) for 
2.1 K. The experimental data is consistent with the ex- 
pected crossover between a lower classical value of (5 at 
high temperature {(5 :^ 0.348 for a 3d-XY universality 
class) and the higher mean field value [5 = 1/2 expected 
for the QPT (T = 0). Note that the effective dimen- 
sionality D = d ^ z > A ioi the spatial dimensionality 
d = 3 and the dynamical exponent z > 1 (2; = 2 for a 
BEC-QCP). 

The magnetic Bragg scattering was measured at three 
different reflections Q = (0.5,0.5,/) with / = 0,3,6 at 
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TABLE I: Experimental magnetic Bragg intensities measured 
at T = 0.16 K and H = 13.5 T and corrected for the Lorentz 
factor, (sin 2^). The calculated intensities for two different 
magnetic structures, one with AI± ||[001] and the other with 
M± II [110] are also listed. 



(0.5, 0.5, 0) 0±2.8 
(0.5, 0.5, 3) 1646.8±32.2 1646.8 
(0.5, 0.5, 6) 55.5±3.6 81.8 



/calc (M^ II [001]) /calc (M^ || [110]) 





1646.8 
327.0 



T = 0.16 K and = 13.5 T. As shown in Table I, their 
relative intensities could be reproduced best by a mag- 
netic structure where the transverse components order 
along the c-axis in the collinear structure that is illus- 
trated in Fig. [21(d). Note that M± is ordered ferromag- 
netically along the [110] direction and antiferromagnet- 
ically along the [100] /[GIG] directions, which is consis- 
tent with the spatially anisotropic exchange interactions 
previously determined by the dispersion relations of the 
triplet excitations with = T [ll[. As shown in the 



upper left inset of Fig. [2](a), for H > Hd the triple ex- 
citation with Sz = I condenses and becomes a gapless 
Goldstone mode with a spin wave velocity of 2.61(14) 
meV/A~^. No anisotropy gap could not be detected 
within the experimental energy resolution of AE = 0.15 
meV, setting the upper limit of the gap size to 0.07 meV. 

In order to find out what kind of perturbations may 
exist in the spin Hamiltonian of this system, we have per- 
formed electron spin resonance (ESR) measurements as 
a function of H and the frequency. These measurements 
probe AS = and A^^ = 1 transition processes [l3 |, 
such as between the spin-triplet states when they are 
split by the external magnetic field. Firstly, in order to 
estimate the gyromagnetic ratio, for Ba3Cr208, the 
ESR measurements were done with a standard DPPH 
(diphenyl-picri-hydrazyl) sample that has the ^-factor 
of 2.00. The ESR transmission spectra are shown as 
solid lines in Fig. ^a). The sharp peak at = 6.8 T 
comes from the standard DPPH while the broad peaks 
at ~ 7 T come from Ba3Cr208. The fact that the 
peak positions of the broad peaks appear at H close to 
but higher than the peak position of DPPH, i^DPPH, indi- 
cates that the ^-factor of Ba3Cr208 is slightly lower than 
2: gab = 2HBPPu/Hab = 1.94 and gc = 2Hbppu/Hc = 
1.93. The order of magnitude of the symmetric exchange 
anisotropy can be estimated as Jo{Ag/g)'^ :^ 1/ieV. 

More ESR measurements with H _L c were done with- 
out DPPH using various frequencies some of which are 
shown in Fig. [3](b). Three ESR modes, labeled by 1, 1', 
and 2, were identified. The mode 2 is due to the ex- 
pected A^* = transitions among the triplet states split 
under H as shown in Fig. El^c), which energy increases 
linearly with H (the black symbols in Fig.[3l^d)). On the 
other hand, for the other two modes, 1 and 1', their char- 
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FIG. 3: (color online) Electron spin resonance transmission 
spectra as a function of external magnetic field (a) at 4.2 K 
with H ± c (blue line) and H \\ c (red line), obtained with 
the frequency of = 190 GHz, and (b) at 1.6 K with H J_ c 
and with two different frequancies u = 149.6 GHz (dotted 
line) and 273.9 GHz (solid line), (c) Energy level scheme 
of the nominal singlet and triplet states as a function of H. 
(d) Frequency dependence of the peak positions of the three 
modes, 1, l\ 2, in H, which are described in the text. 



acteristic energies decrease with increasing following 
hu = huo - giJiBH with vq = 546(9) GHz = 2.2 meV 
(mode 1) and 406(4) GHz = 1.62 meV (mode l^. The i/qs 
correspond to the energies of the singlet-to-triplet exci- 
tations that are 2.22(16) meV and 1.76(17) meV, respec- 
tively, determined by hu;{Q) :^ \/ Jq -\- Jo liQ) where 
7(Q = 0) = 2(J{ + J1 + J'-;') cos(|n^) + 2( + + 
J2 ) + 2(J^ + J'l + Jf )cos(|n7r) with n = and ±1 
for the acoustic and optical modes, respectively, and Js 
are determined by the previous neutron study. These 
indicate that the modes 1 and V are due to the tran- 
sition from the ground state to the | S^Sz >=| 1,1 > 
state. Thus, the ground state of Ba3Cr208 must be a 
mixed state of the singlet and the 6'^ = triplet state, 
I ^ >GS=\ 0,0 > +/ I 1,0 > where / is a mix- 
ing coefficient. Such a mixing is possible if there is a 
Dzyaloshinsky-Moriya (DM) interaction term, D-{siX sj) 
with D \\ H ± c. When H \\ c, the DM term, D ± c 
would mix | 0,0 > with | 1,±1 >, which might be 
the reason for the weak non-zero bulk magnetization for 
H < Hci when H \\ c as shown in Fig. 1 (a). The value 
of M{H - Hci)/M{H = Hc2) < 0.05 should be order of 
I / p-l D/hu{Hci) P where hv{Hci) - 0.4 meV is the 
energy of the ESR signal at yielding | 5 |< 0.1 meV. 
The existence of the DM interactions suggests that the 
inversion symmetry is broken at the center of the dimer 
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FIG. 4: (color online) H — T phase diagram (a) for i7 _L c 
determined from the elastic neutron scattering data (circles) 
shown in Fig. 2 (b) and specific heat data (triangles), and 
(b) for if II c determined from the heat capacity data shown 
in Fig. 1 (b). The lines are fits to a power law (see the text 
for details). 



and the crystal symmetry is lower than the previously 
reported C2/c symmetry 11, iBj . 

The H — T phase diagram was constructed from the 
elastic neutron and specific heat data for if ± c (Fig. [H 

(a) ) and the specific heat data for if || c (Fig. [H 

(b) ). Note that there is only one antiferromagnetic 
phase in both cases. This contrasts to the related com- 
pound Ba3Mn2 08 where, for ii" ± c, a second antifer- 
romangetic phase exists over a narrow region of field, 
which was accounted for by the existence of a single-ion 
anisotropy [l^]. The phase boundary near Hd was fitted 
to a power law, (x| Hc{T) - Hc{0) |^/^. The opti- 
mal fits were obtained with Hci{0) = 11.95(1) T and 
(j) = 2.22(5) for H ± c, and with ifci(O) = 12.13(1) T 
and = 2.21(6) for H \\ c (red lines in Fig. g]). On 
the other hand, when (f) is fixed to the BEG value of 
d/2 = 1.5, the phase boundaries were equally well repro- 
duced up to - 1 K with ifci(O) = 11.92(1) T for if ± c 
and ifci(O) = 12.08(1) T for H \\ c, as shown by the 
blue lines. A previous Monte Garlo simulation study has 
shown that the BEG regime extends up to ~ 0.4 T^ax 
where T^ax is the maximum temperature of the canted 
AFM regime jlO]. This is consistent with our results since 

Tmax ~ 3 K. 

The mixing of the singlet and triplet states in the 
ground state has also been found in other couple dimer 
systems such as KGuGla and SrGu2(B03)2 and 
were ascribed for the existence of DM interactions. The 
DM interaction breaks the U(l) symmetry for if || c and 
consequently it is a relevant perturbation for the BEG- 
QGP. The Ising-like QGP expected for this field orienta- 



tion has an exponent = 2 that should appear below a 
characteristic temperature Tq that depends on the mag- 
nitude of the DM anisotropy. If If Tq <C OATmax one 
should still observe the BEG exponent in the window 
^0 ^ T < OATmax 0- Our analysis of the measured 
Tc{H) curve shows that this data alone is not enough to 
distinguish between the two universality classes. 

In summary, using several different techniques such 
as ESR, specific heat and neutron scattering, we have 
characterized in detail the magnetic interactions and con- 
structed the if — T phase diagram of the newly discovered 
quantum dimer system Ba3Gr208. Our results show that 
near the lower critical field, ifd, the if — T phase dia- 
gram between the disordered and canted AFM states can 
be reproduced by a power law Tc cx| Hc{T) — ifc(O) \^^^ 
that is consistent with (/) = (i/2uptolK~0.4 Tmax- 
Further measurements on other physical quantities near 
the QGP may be needed to single out the universality 
class that dominates in this temperature range. 
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